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I — I Abstract 

K*" We determine the prospects for direct and indirect detection of thermal rehc neutrahnos in 

ly^ supersymmetric theories with muhi-TeV squarks and sleptons. We consider the concrete example 

^-H of the focus point region of minimal supergravity, but our results are generically valid for all models 

with decoupled scalars and mixed Bino-Higgsino or Higgsino-like dark matter. We determine 

^ the parameter space consistent with a 125 GeV Higgs boson including 3-loop corrections in the 

ff^ calculation of the Higgs mass. These corrections increase rra^ by 1-3 GeV, lowering the preferred 

^^ scalar mass scale and decreasing the fine-tuning measure in these scenarios. We then systematically 

K^ examine prospects for dark matter direct and indirect detection. Direct detection constraints do 

S^ not exclude these models, especially for /i < 0. At the same time, the scenario generically predicts 

H spin-independent signals just beyond current bounds. We also consider indirect detection with 

neutrinos, gamma rays, anti-protons, and anti-deuterons. Current IceCube neutrino constraints 

are competitive with direct detection, implying bright prospects for complementary searches with 

both direct and indirect detection. 
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I. INTRODUCTION 

There are now many experimental constraints on weak-scale supersymmetry. These ex- 
clude generic supersymmetric theories in which all superpartners have masses below a TeV, 
and focus attention on the remaining supersymmetric theories that are both phenomenologi- 
cally viable and natural. In this work, we consider focus point supersymmetry [HE], in which 
multi-TeV squarks and sleptons are hierarchically heavier than the other superpartners. 

Focus point models are motivated by a variety of considerations. Heavy first and second 
generation sfermions help satisfy low-energy constraints on flavor and CP violation, and 
heavy third generation sfermions raise the Higgs boson mass to the required level of 125 
GeV [21 Hj. There are also theoretical reasons for expecting scalar superpartners to be 
heavier than the gauginos. For example, such a hierarchy results from an approximate 
U(l)ij symmetry |2j or if none of the supersymmetry-breaking fields is a complete gauge 
singlet [HI E]- Note also that gaugino masses enter the scalar mass renormalization group 
(RG) equations, but scalar masses do not enter the gaugino mass RG equations; letting M1/2 
and mo denote generic gaugino and scalar masses, respectively, the hierarchy ttiq ^ M1/2 
is therefore stable under RG evolution, whereas M1/2 ^ mo is not. Last, although large 
supersymmetry-breaking parameters are generically associated with significant fine-tuning 
of the Higgs potential, simple correlations in high-scale scalar mass parameters may reduce 
the sensitivity of the weak scale to variations in these parameters, providing a naturalness 
motivation for such models. 

In this work, we consider in detail prospects for dark matter detection in such theories [H 
[8]. For concreteness, we consider the focus point region of minimal supergravity (mSUGRA), 
but the results are far more general: when the scalar superpartners are very heavy, they 
effectively decouple from dark matter phenomenology, and the details of the multi-TeV 
spectrum are largely irrelevant. The phenomenology of focus point dark matter encompasses 
the phenomenology of mixed Bino-Higgsino and pure Higgsino neutralino dark matter, and 
our conclusions for dark matter detection are generically valid for any model where the Bino 
soft-supersymmetry breaking mass is lower than the Wino mass. 

In Sec. [TT| we explain our treatment of mSUGRA parameter space. We then turn to the 



Higgs mass in Sec. HI, including a 3-loop calculation of the Higgs mass using the public 
code H3m [9l [10]. We find that 3-loop contributions raise the Higgs mass by 1-3 GeV over 
2-loop results. Given the logarithmic sensitivity of the Higgs mass to the top squark mass, 
this lowers the preferred range of stop masses considerably. In this calculation stop masses 
as low as 3 TeV are consistent with the measured Higgs mass, even without significant stop 
left-right mixing. In the focus point parameter space, this correlates with a gluino as light 
as 2 TeV. 

We then consider prospects for dark matter detection in the region of parameter space pre- 
ferred by the Higgs mass and other phenomenological constraints, including direct searches 
for supersymmetric particles. In Sec. [TV] we discuss both spin-independent and spin- 
dependent direct detection and show that, contrary to claims in the literature, perfectly 
viable regions of parameter space remain, especially for /i < 0. Crucial to this conclusion is 
the small value for the strange quark content of the nucleon now preferred by both lattice 
calculations and chiral perturbation theory results. At the same time, the scenario generi- 
cally predicts spin-independent cross sections a?^ ~ 1 zb = 10~^ pb = 10~^^ cm^, implying 
that dark matter candidates in this class of theories might very well be discovered by direct 
detection experiments in the near future. 



In Secs.jV, VI, and VII we analyze the implications for indirect detection with neutrinos, 



gamma rays, and anti-matter, respectively. Although gamma rays and anti-matter are 
currently not very constraining in focus point scenarios, current bounds from observations 
of neutrinos from the direction of the Sun with IceCube are stringent, and future runs with 
planned upgrades will probe much of the preferred region, providing an exciting, and in 



many respects orthogonal, complement to direct detection. In Sec. |VIII[ we discuss our 
results and conclude. 



II. PARAMETER SPACE AND LHC SUPERPARTNER SEARCHES 

The defining feature of focus point supersymmetry is the insensitivity of the weak scale 
to variations in the fundamental supersymmetry-breaking parameters, even in the presence 
of multi-TeV soft supersymmetry-breaking parameters. Focus point supersymmetry accom- 
modates a range of thermal relic neutralinos that vary continuously from ~ 100 GeV Bino- 
Higgsino mixtures to heavier and more Higgsino-like neutralinos, culminating in Higgsino- 
like neutralinos with masses around 1 TeV [7J . Given the appeal of neutralino dark matter, 
it is natural to impose the thermal relic density as a constraint on the parameters space. In 
the context of mSUGRA, this constraint allows for a departure from the typical (mo, M1/2) 
parameter space — in which the cosmologically viable region is only a small sliver — to a 
parameter space in which every point is cosmologically viable and more parameters can be 
examined [TT]. This parameter space is particularly relevant in light of the first three years 
of LHC results, which have effectively eliminated the so-called "bulk" scenario for neutralino 
dark matter with light scalars and severely constrained coannihilation scenarios with light 
scalars, while leaving the focus point relatively unscathed and strong as a possibility for 
neutralino dark matter. 

In mSUGRA, the relic density constraint can be cast as the requirement that 

Q^ (mo, Ml/2, ^0, tan (3, sign(/i)) = ^dm , (1) 

where ^dm — 0.23 is the dark matter density in units of the critical density. Focus point su- 
persymmetry is possible with large A-parameters [12], but given the motivations of simplicity 
and a hierarchy between supersymmetry-breaking parameters enforced by an approximate 
U{1)r symmetry, we choose Aq = throughout. We may then use Eq. ([I]) to solve for mo 
and present results in the (tan/3, M1/2) plane for both signs of fi, with every point in these 
planes having the correct relic density. In general, Eq. ([I| may be satisfied by more than 
one value of niQ; for example, there may be a coannihilation solution at low ruQ and a focus 
point solution at larger mo. In such cases, we always use the largest allowed value of mo. 

Figure [l| shows contours of mo in the (tan/3,Mi/2) for points satisfying the relic den- 
sity constraint, using SOFTSUSY 3.1.7 [14] to generate the SUSY spectrum and Mi- 
crOMEGAs 2.4 [^ to calculate the relic density. These solutions for mo are found for 
low values of |/i| located near the //^ < region, where radiative electroweak symmetry 
breaking fails. The /i^ < region moves to higher mo ior increasing M1/2 and decreasing 
tan P due to RG effects, and this behavior is reflected in the mo contours. In Fig. [T] the 
shaded region with low M1/2 is excluded by ATLAS searches for jets + missing energy [T3] . 
The other shaded regions, which will appear in all of our figures, include a region at large 
tan/3, where the RG evolution in SOFTSUSY becomes unreliable, and a region at large 
Ml/2 for n < 0, where numerical issues with loop corrections to neutralino masses make 
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FIG. 1. Contours of mo in TeV (solid blue) and fine-tuning parameter c (dot-dashed gold) in the 
(tan/3,Mi/2) plane for Q^ ~ 0.23, Aq = and /i < (left) and // > (right). The red shaded 
regions at low M^/2 are excluded by the ATLAS gluino bound [13]. In the gray shaded regions at 
large tan/3, the RG evolution in SOFTSUSY becomes unreliable, and in the green shaded region 
at large M1/2 for /i < 0, numerical issues with loop corrections to neutralino masses make the 
solution algorithm for il unreliable. 



the solution algorithm for ft unreliable. We stress that these last two regions are excluded 
not by theoretical or experimental constraints, but rather because numerical complications 
hinder our ability to make accurate predictions. 

Since the sfermion sector is decoupled in focus point supersymmetry, the properties of 
neutralino dark matter are determined primarily by its mass and the amount of Bino- 
Higgsino mixing present. If the gauge eigenstate composition of the lightest neutralino is 
given by 
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with a-Of -C 1 in the focus point region, the dominant processes for both annihilation and 



scattering are proportional to either 
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mixing can be usefully parameterized by the Bino content a^. Figure [2] contains contours of 
m^ and a^ consistent with VL^ = Qbm- For much of the parameter space, the neutralino dark 
matter is a Bino-Higgsino mixture, but as M1/2 increases, m^ increases, and a^ decreases: 
the increasing Higgsino content compensates for the suppression of the annihilation cross- 
section by larger neutralino masses to keep the thermal relic density constant. The behavior 
is similar for both signs of /x, though a^ is somewhat larger in the /z < case relative to the 
/i > due to the relative signs of a^ for different signs of /i. In the limit of large M1/2, 
the neutralino becomes nearly pure Higgsino with a^ — )■ 0, and the neutralino mass reaches 
rriy. ?a 1 TeV. 

In focus point scenarios, the weak scale is relatively insensitive to variations in supersym- 
metry breaking parameters, allowing for improved naturalness even with multi-TeV sfermion 
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FIG. 2. Contours of m^ in GeV (black dotted) and |a^| (solid colored). 

masses. There are many prescriptions for quantifying this naturalness, all of which are sub- 
ject to significant subjective choices; for a review, see Ref. [ISj. Here we use a naturalness 
measure based on the sensitivity coefficients [T71 ITS] 
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where a^ is one of the input GUT-scale parameters ml, M^^? ^O' /^O' ^^^ "^3' ^^^ ^u^d 
soft mass parameter. The overall fine-tuning of a model is defined as 



c = max{ca} 



(4) 



and contours of c are shown in Fig. [T] In the explored region, Cmo is always the largest 
sensitivity coefficient, and contours of c roughly follow contours of mo, with values of m,o ~ 
4 TeV corresponding roughly to c ~ 250. A subset of the mSUGRA boundary conditions 
implies focusing, and the values of c shown in Fig. [T] are much smaller than would be expected 
without the focus point behavior. 

III. CONSTRAINTS FROM THE HIGGS MASS 



The mass of the recently-discovered SM-like Higgs boson P, \^ provides a stringent con- 
straint on the parameter space of any supersymmetric model. The most recent mass mea- 
surements are [191 120] 



ATLAS M : 124.3 +°J +°j GeV 
ATLAS 77 :126.8 ± 0.2 ± 0.7 GeV 
CMS 4i :125.8 ± 0.5 ± 0.2 GeV 
CMS 77 :125.4 ± 0.5 ± 0.6 GeV 



(5) 
(6) 
(7) 
(8) 




1 .5 -- 



0.5 



129 

2o Theory Error 
■,;,-:::::::::::::=:j^--^^- 
lo Theory Error 
1^5 



121 



123 



:dn 



dj 



dj 



■,119, 



*; 
o 

CO 



o 

CL 



10 15 20 25 30 35 40 45 50 55 60 

tanp 

(a) /z < 



2.5 



^ 



1.5 



0.5 



- 




ri 


- 




i 


""'•.,_ 




3. 


- 


129 


>; 






a 


- 


2a Theory Error 


o 

If) 














. ' 


i27 


c 






Q. 


- 


la Theory Error 


T3 
1 


- 


125 




123 


,1,,,, 




121 


.(.r-rM't"i"i"'i I 


1 [ [ : |1T9'7 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 ; 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



10 15 20 25 30 35 40 45 50 55 60 

tanp 

(b) /x > 



FIG. 3. Contours of rrih in GeV. In the shaded regions, the theoretical prediction for rri}^ is within 
the indicated theoretical uncertainty of the experimental central value ruh = 125.5 GeV. 



where the first uncertainties are statistical and the second uncertainties are systematic. 

We calculate the lightest Higgs mass in the focus point region of mSUGRA with the pro- 
gram H3m, which calculates nth in the !5r scheme including the dominant 3-loop contribu- 
tions at 0{ata1) [9], [10]. In addition, we modified H3m to increase the precision in the calcu- 
lation of the running DR top quark massjj We set m^° "^ = 173.2 GeV and as{mz) = 0.1184, 
and fix the renormalization scale to the geometric mean of the stop masses. For further 
details, see Ref. pi] . 

In Fig. |3} we plot contours of rrih in the parameter space defined by Fig. [ij We find that 
the 3-loop terms generate a 1-3 GeV increase in rrih over the 2-loop truncation. The 2-loop 
terms in turn generate a 5-8 GeV increase over the 1-loop truncation, indicating convergence 
of the series. We observe also that the improved treatment of mf^ and a^^ increases the 
2-loop prediction relative to FeynHiggs. For comparison, note that the geometric mean of 
the stop masses ranges from about 1 TeV at low M1/2 to 8 TeV at high M1/2 in the plotted 
parameter space. 

In Fig. [3] we shade regions where the difference between the calculated rrih and the tenta- 
tive central value 125.5 GeV is within the indicated theoretical uncertainty in the calculation. 
At each point on the parameter space, we assign a theoretical error bar Ath, defined as 



A. 



th 



^pert — 



(Apcrt) + (Apara) 



m 



(3-loop) 
h 



m 



(2-loop) 
h 



A 



para 



mh{mt = 174.2 GeV) - mh{mt = 173.2 GeV) . 



(9) 



^ These changes are incorporated in the next version of H3m, which will be released simultaneously with 
Ref. [21]. 



The uncertainty Apcrt from higher-order terms in the perturbation series is estimated to be 
in the range 0.5-1.5 GeVnThe parametric uncertainty Apara induced by the uncertainty in 
the top quark mass is typically of order 0.5-1 GeV in the focus point parameter space. 

The positive 3-loop terms significantly impact the preferred range of superpartner masses. 
Requiring that the theoretical prediction be within Ath of the experimental result and 
neglecting the significant experimental uncertainties, scalar mass parameters as low as 
mo ~ 4 TeV, corresponding to stop masses as low as 3 TeV, and gluino masses as low 
as rrig ~ 2.8M1/2 ~ 2 TeV are consistent with the measured Higgs mass. Note that, combin- 
ing the results shown in Figs. [T] and [3} the 3-loop rrih contributions also significantly decrease 
allowed values of the fine-tuning parameter c. 

IV. DIRECT DARK MATTER DETECTION 

It is well-known that thermally-produced neutralinos can possess a wide range of direct 
detection cross sections, from those that are significantly excluded to those that are orders 
of magnitude below current sensitivities. However, this full range of cross sections is not 
generic. Highly suppressed direct detection is typically associated with pure Bino scenarios, 
which have the correct thermal relic density only if there are light sfermions, co-annihilation, 
or resonant annihilation through the pseudoscalar Higgs resonance. The first two possibilities 
are disfavored by the non-observation of light squarks at the LHC, while the third depends 
upon careful tuning of the pseudoscalar Higgs mass to ttia ^ 2m^. Most of the remaining 
parameter space is populated by models with Bino-Higgsino mixing like that found in the fo- 
cus point region. For these models, the Bino-Higgsino mixing also sets the spin-independent 
neutralino-proton scattering cross section, which falls in the range a?^ ~ 1 — 40 zb for a wide 
range of model parameters when neutralinos have the right thermal relic density [H] . This 
range of cross sections is particularly relevant for current and near-future direct detection 
experiments; the XENONIOO experiment [231 El] has begun probing this range of relevant 
cross sections, and near-future direct detection experiments will be sensitive to most of the 
focus point region of mSUGRA. 

In the focus point region, (7^^ is dominated by Higgs-mediated diagrams, and the Higgs- 
neutralino coupling is sensitively dependent on the sign of /i, producing a suppression of 0"^^ 
in the /i < case relative to the /x > case. For moderate tan/3 this leads to a relative 
factor of a few in a?^, from the coupling coefficients and at large tan/3 due to the relative 
contribution of the heavy Higgs-mediated diagrams. Although the general lore holds that 
yU > is preferable to address the discrepancy in the anomalous magnetic moment of the 
muon [25H27] . in focus point theories the contribution for either sign of /i is too small to 
produce consistency or to further aggravate the discrepancy without considering significant 
non-universality of smuon masses [28] . 

Determinations of a?^ for neutralinos also suffer from the well-known uncertainty in the 
quark scalar form factor of the nucleons, f^ , defined as 

{N\m,tPgtlj,\N)=f^MN . (10) 

The form factors for the up- and down-type quarks are well measured, and the heavy quark 
form factors are determined by loop contributions from the gluon form factor, but there is a 



^ The size of the 3-loop corrections is consistent within the uncertainty with the NLL analysis of Ref. [22] 
which used a somewhat different organization of the perturbation series. 




10 15 20 25 30 35 40 45 50 

tanp 

(a) /i < 



55 60 




15 20 25 30 35 40 45 50 55 60 



tanp 



(b) /x > 



FIG. 4. Contours of spin-independent scattering cross section 0"^^ in zb. The shaded regions are 
excluded by XENONIOO [Mj, assuming local dark matter density piocai = 0.3 GeV/cm^ (light 
shaded) and piocai = 0.15 GeV/cm^ (dark shaded). 



longstanding controversy regarding the strange quark form factor, which feeds into 0"^ in a 



quantitatively important way [29143T] . Older results from chiral perturbation theory 
combined with determination of the nucleon sigma term from meson scattering data 
and supported by direct computation [3S], suggested fs = /" = /f ~ 0.36. For this value 
of fs, the other form factors are all much smaller, /j^^ < 0.05, and so the strange quark 
contribution dominates the direct detection cross section [29]. However, recent lattice studies 
favor a much smaller value of [371 EHI 



/. ~ 0.05 



(11) 



much closer to the other quark flavors [30', 13^. It has also been argued that the lower value 
for fs is consistent with chiral perturbation theory computations, provided higher-order 
baryon decuplet contributions are taken into account [36l|37l HOl HI] . A recent calculation 
considering these contributions found fs = 0.017 ± 0.15 [12]; for similar recent conclusions, 
see Refs. [121 HI]. Here we take fs = 0.05 in deriving direct detection cross sections. 

Figure [4| shows exclusion contours for XENONIOO in the (tan/3, M1/2) plane for both 
signs of /i and piocai = 0.3 GeV/cm^. For /i > 0, current XENONIOO bounds require 
^1/2 ^ 1-8 TeV for a wide range of tan/3, with stronger exclusions at large and small tan/3, 
as discussed above. A small region at very large tan/3 is allowed for M1/2 ^ 500 GeV; 
here the lightest neutralino is nearly pure Bino due to the A-funnel crossing through the 
focus point region. For /i < 0, XENONIOO requires M1/2 > 1.3 TeV for moderate values 
of tan/3, but the exclusions are much weaker for small and large tan/3, and large portions 
of the parameter space remain viable. Exclusion contours for piocai = 0.15 GeV/cm^ are 
also presented, motivated by the possibility of a local dark matter density somewhat lower 
than normal due to the presence of small-scale structure [45j. For this lower value of piocai 
and both signs of /x, the excluded region is roughly comparable to that excluded by gluino 
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FIG. 5. Contours of the spin-dependent neutralino-proton scattering cross section Up in pb. The 
shaded region indicates the reach of COUPP-60 [l6] after 12 months with piocai = 0.3 GeV/cm^. 



searches, only becoming stronger for large and small tan/3 when /x > 0, and almost none of 
the parameter space preferred by the Higgs mass is excluded by direct detection. 

Dark matter may also be detected directly through its spin-dependent couplings. Fig- 
ure p^ shows contours of constant 0"^°, the spin-dependent neutralino-proton scattering cross 
section. Across the parameter space of the focus point region compatible with the correct 
thermal neutralino relic density, a?^ is in the range 10~^ — 3 x 10~^ pb, for both signs of 
fi, decreasing with increasing M1/2. At large values of M1/2, the lightest neutralino becomes 
increasingly Higgsino-like, suppressing cjp^. However, the observed Higgs mass disfavors the 
pure Higgsino limit, and the 2a allowed region for rrih favors cr^^ in the range 10"'' — 10~^ 
pb. 

The shaded region shows the sensitivity expected from COUPP-60 [IHl ST] , corresponding 
to a data-taking period of 12 months at SNOLAB, in the zero-background assumption and 
using the typical local dark matter density of piocai = 0.3 GeV/cm^. With one year of data, 
the COUPP-500kg experimental sensitivity is anticipated to range between a few x 10^^ pb 
at 100 GeV to a few x 10~^ pb at 1 TeV, thus covering a significant portion of the parameter 
space of interest here. 



V. NEUTRINOS FROM ANNIHILATION IN THE SUN AND IN THE EARTH 



The search for high-energy neutrinos from the direction of the center of the Sun or of the 
Earth has a special place in the ranks of indirect detection techniques. In the limit where the 
capture rate of dark matter particles in celestial bodies is equilibrated by the annihilation 
rate, the flux of neutrinos solely depends on the scattering cross section of dark matter off 
of nuclei in the celestial bodies. In the case of the Sun, the dominant scattering mechanism 
for neutralinos in the minimal supersymmetric standard model is typically spin- dependent 
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FIG. 6. The flux of muons in units of km yr ^ with energies above 1 GeV at IceCube. The 
shaded region is excluded by current limits from IceCube/DeepCore 



scattering, while scattering in the Earth is dominated by spin-independent processes. Unlike 
searches for antimatter or gamma rays, where the target dark matter densities are generally 
poorly known and affected by large uncertainties, the flux of neutrinos from the Sun or the 
Earth has a rather mild dependence on astrophysical inputs. The only crucial information 
is, in fact, the local dark matter density. In this respect, of all indirect searches, neutrino 
telescopes provide perhaps the most robust limits. 

In Fig. |6] we show the flux of muons produced via charge-current interactions by high- 
energy neutrinos from dark matter annihilation in the Sun. To calculate this rate (as well 
as all of the subsequent indirect detection rates) we employ the DarkSUSY package, version 
5.0.5 |l9]. Figure [6] shows the integrated muon rate for muons with energies larger than 1 
GeV. The shaded region at the bottom is excluded by the latest results from 317 days of data 
taken from 2010-11 at the IceCube neutrino telescope with the 79-string conflguration, and 
with the use of the DeepCore sub-array [48j . This region excludes a parameter space portion 
comparable to that excluded by current LHC searches. Note that the 1 GeV threshold is 
much lower than the detector's actual energy threshold, even with the use of DeepCore, 
but the 1 GeV threshold is used in Ref. |18] for consistency with other results in the fleld, 
especially from experiments such as SuperKamiokande, where the 1 GeV threshold is actually 
experimentally meaningful. For IceCube/DeepCore, the extrapolation below the 1 GeV 
threshold is made based on the assumed neutrino spectrum, which in the focus point region 
corresponds closely to the W^W~ channel for which the exclusion limits are quoted in 

Ref. m- 

The results of Ref. [IE] fell short by about a factor 2-5 of the anticipated target sensitivity 
quoted in Ref. [50] for 180 days. We flnd that had the detector performed to the level 
anticipated in Ref. [50], the exclusion limit would have extended up to M1/2 ~ 1-5 TeV, 
covering much of the parameter space of the focus point region compatible with the Higgs 
mass. This is supported by Fig. [7], where we show the flux of muons integrated above 
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FIG. 7. The flux of inuons in units of km" yr~^ with energies above 100 GeV at IceCube. 
The shaded region shows the originaUy anticipated sensitivity region for 180 hve days for the 
IceCube/DeepCore system [50] , 



a 100 GeV threshold; these numbers are therefore more indicative of the actual number 
of events IceCube might detect than those shown by the contours of Fig. [6] The shaded 
region corresponds to the original 180 days sensitivity target, which would have excluded 
Ml/2 ^ 1-5 TeV with little dependence on tan/3, corresponding to a lightest neutralino mass 
of ~ 600 GeV. This emphasizes how promising neutrino telescope searches are in the context 
of searches for a signal of new physics from the focus point region. We also note that the 
recent null results from the ANTARES collaboration [5Tj reinforce the lack of a high-energy 
neutrino signal from the Sun, at a level very close to the current IceCube/DeepCore limits. 
The rates of high-energy neutrinos, and consequently of muons, from neutralino annihi- 
lation in the center of the Earth are not nearly as exciting as those from the center of the 
Sun. We show in Fig. |8] the calculated fluxes of muons from the Earth, again integrated 
above a 1 GeV energy threshold. Nowhere do we obtain fluxes much larger than 10~^ km~^ 
yr~^, which is clearly well below the sensitivity of km^-sized neutrino telescopes. We note 
that unlike the case of the Sun, for the Earth the dependence of the flux of neutrinos on the 
spin-independent cross section induces a signiflcant dependence on the sign of /i, with pos- 
itive n producing larger fluxes due to the lack of interfering terms in the neutralino-proton 
scalar cross section, as discussed in the previous section. 



VI. GAMMA RAYS 



Gamma rays provide another promising possibility for the indirect detection of dark 
matter. This signal is especially relevant now that the Fermi Large Area Telescope (LAT) [52] 
has revolutionized our understanding of the high-energy sky, in a photon energy range 
extraordinarily relevant for indirect searches for WIMP dark matter. 
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FIG. 8. The flux of muons in units of km yr ^ with energies above 1 GeV at IceCube from dark 
matter annihilation in the center of the Earth. 

The gamma-ray signal may take one of two forms. It may appear as a monochromatic 
line, if photons are produced as one or both of the annihilation products in a two-body final 
state. Alternatively, the signal may be an excess of continuum gamma rays extending for 
several decades in energy below the dark matter particle mass. Such continuum gamma 
rays are typically produced from the two-photon decay of neutral pions resulting from the 
hadronization of annihilation products, or from final state radiation, or from inverse Comp- 
ton processes associated with final state electrons and positrons. 

We begin by considering the line signal. Figure [9] shows curves of constant branching ratio 
into two photons. The branching fraction increases towards increasing massesjjbut is always 
much smaller than the percent level. In the focus point parameter space, the thermally- 
averaged neutralino pair annihilation cross section always lies at about 1 — 2 x 10^^^ cm^ s~^, 
with the mismatch with the canonical value of 3 x 10^^^ cm'^ s~^ being due to chargino 
and next-to-lightest neutralino coannihilation. These values imply that the Fermi LAT 
Collaboration line search limits [54j do not yet constrain this parameter space. 

The recent discovery of a 130 GeV line-like feature in the Fermi LAT data has attracted 
great attention [551 ES]- Our results indicate that focus point supersymmetry does not 
provide a viable framework to explain the line feature with dark matter annihilation, as 
the branching ratio into two photons, and the associated pair-annihilation cross section, are 
much smaller than the required value of ~ 10~^^ cm^ s~^. 

Turning next to the continuum signals, we consider annihilation in local dwarf galaxies, 
currently one of the most stringent and robust limits on the pair-annihilation cross section 
of dark matter. Cross sections of the order of what the theory predicts over the parameter 
space of interest are only constrained for neutralino masses on the order of 30 GeV [57]. In 

^ This might be partly due to the fact that the annihilation mode into two photons is the only electroweak 
one-loop correction implemented in DarkSUSY: this artificially boosts the branching ratio into two photons 
as the neutralino mass approaches M^/aw jS'Sj : this result should therefore be taken with a grain of salt. 
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FIG. 9. Gamma ray searches in the focus point region. The curves indicate constant values for the 
branching ratio for neutrahno pair annihilation into two photons. Null results from Fermi searches 
for a monochromatic gamma-ray line do not put any constraint on this plane. Null results for 
continuum gamma-ray signals with Fermi using stacked dwarf galaxies also do not exclude any of 
this parameter space. However, improvements of current bounds on the gamma-ray continuum 
will probe the parameter space. The shaded regions indicate the performance of searches for a 
continuum gamma-ray signal, assuming current sensitivities are improved by factors of 5 and 20, 
as indicated. 



focus point supersymmetry, such masses are never consistent with the relic density constraint 
(and are also excluded by neutrino telescope searches and by LHC results), and the limits 
weaken approximately quadratically with mass. 

This is illustrated with the shaded regions shown in Fig. |9j which indicate the improve- 
ment to the Fermi limits needed to probe the parameter space of interest; we indicate the 
sensitivity lines corresponding to improvements by factors of 5 and 20. To approach the 
level of Ml/2 ~ 1-5 TeV, the Fermi limit from stacked dwarf galaxies [SZ] would need to be 
improved by a factor of 20. Such an improvement would take a time- frame which is beyond 
the anticipated lifetime of the mission. 

It is important to note that we have considered here line and continuum signals given 
conservative assumptions. Constraints can be obtained by employing optimistic choices for 
the density profile of the inner Galaxy, or of external galaxies or clusters, or by utilizing 
optimistic assuptions for the dark matter sub-structure content and structure. Here, we 
have limited ourselves to the more conservative limits obtained by the Fermi Collaboration 
for line signals [51] and continuum signals from stacked dwarfs [2Z] . We emphasize that had 
we used the Galactic center and a favorable dark matter density profile, we could have easily 
reached radically more optimistic conclusions. 

We do not show here predictions for the performance of a future Cherenkov Telescope 
Array (CTA); see, e.g., Ref. |58]. Certain sensitivity estimates for the reach of CTA optimisti- 
cally carve into the parameter space of the focus point region, for example from observations 
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FIG. 10. The differential flux of anti-protons in units of GeV cm ^ s ^ sr ^ at an energy of 19.6 
GeV. 

of the inner Galaxy [5SJ • Interestingly, CTA will be especially sensitive to WIMP masses in 
the TeV region, and is thus, in principle, an ideal instrument to look for a signal in the focus 
point region. Under conservative assumptions, however, CTA, like Fermi, is not guaranteed 
to detect a signal from dark matter models in the focus point region. In addition, annihila- 
tion of a 1 TeV neutralino in the focus point region to the level needed for a detection with 
CTA would lead to significant low-energy inverse Compton gamma-ray production, which 
might conflict with existing Fermi LAT limits. We postpone detailed discussion to future 
work, but we emphasize that CTA will be a key observational tool in the search for particle 
dark matter in this region, especially if a signal for TeV-mass dark matter were detected in 
direct detection or neutrino telescope experiments. 



VII. ANTIMATTER 



The successful deployment of the Alpha Magnetic Spectrometer (AMS-02) on board the 
International Space Station has boosted hopes and expectations of using cosmic-ray anti- 
matter as a probe of annihilation of Galactic dark matter. In the context of the focus point 
region, for models with the correct thermal neutralino relic density, the flux of positrons is 
always too small to be detectable with any significance by current experiments, so we focus 
here on anti-protons and anti-deuterons. The latter choice is motivated by the extremely 
suppressed background rate and great discrimination capabilities against anti-protons that 
the future General Antiparticle Spectrometer (GAPS) mission promises for anti-deuterons 
in the low energy (approximately at or below 1 GeV) range 



Figure [TO] shows the flux of anti-protons at an energy of 19.6 GeV. We use the default 
propagation parameters for charged cosmic rays in DarkSUSY, as well as the default Dark- 
SUSY |19] dark matter halo density profile. We choose the particular energy of 19.6 GeV 
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FIG. 11. The flux of anti-deuterons at an energy of 1 GeV, in units of GeV cm ^ s ^ sr ^. 



for two reasons: 

(1) It was shown in Ref. ^Glj (see Fig. 10, left and Fig. 11, left) that the best signal to 
background ratio for anti-proton searches in the focus point region ranges between 10 and 
100 GeV in kinetic energy, with an optimal value of about 20 GeV when factoring in the 
need to observe a large enough number of signal events; 

(2) 19.6 GeV corresponds to the central value of the relevant energy bin reported by 
the PAMELA Collaboration ^. At that energy, PAMELA quotes a flux of 7.2 x 
10-s GeV"^ cm-2 s'^ st-\ 



The contours in Fig. 10 indicate that the signal-to-noise ratio expected in the "sweet 
spot" for the anti-proton kinetic energy ranges between 2% for very light neutralinos to less 
than 0.1% for more massive neutralinos. We flnd almost no variation between negative and 
positive values of fi. Given the absence of any striking spectral feature in the predicted 
spectrum of anti-protons in the focus point region [HI], and the fact that variations in the 
cosmic ray anti-proton diffusion and energy loss parameters can induce deformation to the 
background spectrum much larger than the percent level, we conclude that the predicted 
flux of anti-protons is generically too small to provide a conclusive dark matter detection 
avenue for neutralinos in the focus point region. 

Figure [TT] shows the anti-deuteron flux at 1 GeV. Although the GAPS experiment will 



primarily target lower energies (likely between 0.1 and 0.3 GeV), the AMS-02 limits are 
likely to be best in the 1 GeV range. In addition, the predicted flux at 0.1-0.3 GeV is 
typically comparable (within less than a factor 2) to that at 1 GeV for neutralinos in the 
focus point region. 

The GAPS experimental sensitivity target is at present estimated to be at the level of 
just under 10~^^ GeV~^ cm^^ s^^ sr^^, while AMS-02 should be able to reach a sensitivity 
of about 10~^° GeV~^ cm~^ s^^ sr^^, or approximately an order of magnitude less constrain- 



ing than GAPS. Figure [TT] therefore illustrates that across the relevant focus point region 
parameter space the expected anti-deuteron signal is between 1 and 4 orders of magnitude 



15 



smaller than the best foreseeable experimental sensitivity, making this indirect detection 
channel inconclusive to search for a dark matter signal. 



VIII. DISCUSSION AND CONCLUSIONS 

In the context of supersymmetric extensions of the standard model, the discovery of a rel- 
atively heavy Higgs boson at the LHC, coupled with null results from superpartner searches, 
provides strong motivation for considering models with multi-TeV squarks and sleptons. We 
consider cosmologically-motivated focus point model realizations of this scenario in which 
dark matter is entirely composed of thermal relics that are mixed Bino-Higgsino or Higgsino- 
like neutralinos. 

Our main findings are the following: 

• These models remain viable. Claims to the contrary are apparently the result of (a) 
requiring supersymmetry to resolve the {g — 2)^ discrepancy (a requirement that is 
tantamount to considering the standard model to be excluded by this discrepancy), 
(b) considering only /i > (presumably for historical reasons linked to (a)), (c) using 
large values of fs that are now highly disfavored, (d) imposing some highly subjective 
naturalness criterion, or (e) a combination of these. 

• The leading 3- loop 0{ata'^) contributions to the Higgs mass are positive, lowering the 
preferred values of scalar masses (possibly to values within reach of the LHC) and 
improving the fine-tuning of these scenarios. 

• Some focus point parameter space is excluded by bounds from direct searches for dark 
matter, but some remains, including much of the parameter space with n < 0. In the 
allowed regions, the predicted spin-independent cross sections are just beyond current 
bounds from XENON, and spin-dependent scattering is also close to the experimental 
sensitivity expected in the near future. 

• For indirect detection, searches for neutrinos from the core of the Sun at Ice- 
Cube/DeepCore exclude focus point neutralinos lighter than about 170 GeV. The an- 
ticipated detector performance would have placed constraints on neutralinos as heavy 
as 600 GeV, covering most of the focus point parameter space. There are therefore 
bright prospects for dark matter discovery through neutrinos at IceCube/DeepCore. 
Similar sensitivity is being reached by other experiments, such as ANTARES. These 
results are insensitive to halo model choices, and also do not depend on, e.g., the 
strange quark content of the proton, and so yield promising probes that are highly 
complementary to direct detection. 

• The predicted neutrino flux from the center of the Earth is many orders of magnitude 
below detectability. 

• We have also considered gamma rays from the galactic center and dwarf galaxies 
producing either line or continuum signals. Signals in gamma rays are not as promising 
as in neutrinos, at least for the conservative choices of the relevant dark matter density 
proflles we employed here, but may nevertheless still be seen in future experiments such 
as CTA. 
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• For indirect detection of anti-protons, the signal-to-background ratio, even at optimal 
energies, is at the percent level and too small to provide a convincing avenue for dark 
matter detection. 

• Anti-deuteron rates are one to four orders of magnitude below the foreseen experimen- 
tal sensitivity of future dedicated experiments, such as GAPS. 

To summarize, LHC results so far motivate focus point supersymmetry, which has excit- 
ing implications for dark matter searches. Among the most promising are direct searches 
for spin-independent scattering and indirect searches with neutrino telescopes, but other ap- 
proaches discussed here may also yield signals. Uncertainties in the Higgs mass calculation 
also leave open the possibility that squarks and gluinos may be within reach of the LHC, 
even without large left-right stop mixing. If focus point supersymmetry is realized in nature 
and focus point neutralinos make up all of the dark matter in the Universe, a signal in one 
or more of the complementary probes (colliders, direct detection, and indirect detection) 
will appear in the coming few years. 
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